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Abstract 

It has been known for several decades that Einstein's field equations, when 
projected onto a null surface, exhibits a structure very similar to non-relativistic 
Navier-Stokes equation. I show that this result arises quite naturally when 
gravitational dynamics is viewed as an emergent phenomenon. Extremising 
the spacetime entropy density associated with the null surfaces leads to a set 
of equations which, when viewed in the local inertial frame, becomes identical 
to the Navier-Stokes (NS) equation. This is in contrast to the usual description 
of Damour-Navier-Stokes (DNS) equation in a general coordinate system, in 
which there appears a Lie derivative rather than convective derivative. I discuss 
this difference, its importance and why it is more appropriate to view the 
equation in a local inertial frame. The viscous force on fluid, arising from 
the gradient of the viscous stress-tensor, involves the second derivatives of the 
metric and does not vanish in the local inertial frame while the viscous stress- 
tensor itself vanishes so that inertial observers detect no dissipation. We thus 
provide an entropy extremisation principle that leads to the DNS equation, 
which makes the hydrodynamical analogy with gravity completely natural and 
obvious. Several implications of these results are discussed. 



1 Gravity as an emergent phenomenon 

Large amount of theoretical evidence suggests that gravity could be an emergent 
phenomenon like gas dynamics or elasticity with the gravitational field equations 
having the same status as, say, the equations of gas/fluid dynamics P^. They seem 
to describe the thermodynamic limit of the statistical mechanics of (as yet unknown) 
atoms of spacetime or, equivalently, hydrodynamic (long wavelength) limit of a suit- 
ably defined 'fluid' system. 

A very direct argument in favour of this paradigm is the validity of principle of 
equipartition which allows one to use the classical field equations, along with the 
expression for Davies-Unruh temperature [5], to determine the number density of 
microscopic degrees of freedom [5]. One finds that any null surface in the space- 
time is endowed with a number density of microscopic degrees of freedom ('atoms 
of spacetime') which depends on the structure of the gravitational theory. In the 
simplest context of Einstein's theory, this density is a constant (equal to one degree 
of freedom per Planck area) while in a generic Lanczos-Lovelock model, it is propor- 
tional to P'^^^'^ = dL I dRabcd where L is the Lanczos-Lovelock Lagrangian and Rabcd 
is the curvature tensor [3]. We shall hereafter call P'^^^^ as the spacetime entropy 
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tensor or simply the entropy tensor of the theory. (In the case of Einstein's theory 
with L (X R, we have P^^^ oc {S^S'^ — S'^S^)). Given the number density of microscopic 
degrees of freedom one can construct an entropy functional for the system which 
will depend only on P'^'"^'^ and lead correctly to the Wald entropy [1] of the Lanczos- 
Lovelock theory. The crucial fact that entropy of a horizon in a gravitational theory 
depends on the theory (and is proportional to the area only in the simple case of 
Einstein gravity) finds expression in the related fact that the density of microscopic 
degrees of freedom determined through the law of equipartition, depends on the 
theory. 

This key 'internal evidence' from the structure of gravitational theories — es- 
pecially in Lanczos-Lovelock models — suggests that they may be viewed as the 
thermodynamic limit of the statistical mechanics of the microscopic degrees of free- 
dom. This is reminiscent of the work by Boltzmann and others, who used the 'in- 
ternal evidence' of thermal phenomena to deduce the existence of atoms/molecules 
and determine the Avogadro number which quantifies their number density even 
without knowing what they are and having no direct observational support for their 
existence. 

The above results were obtained by starting from the field equations of the theory, 
rewriting them in the form of law of equipartition and thus determining the density 
of microscopic degrees of freedom. A more attractive procedure will be to invert this 
argument and obtain the field equations of the theory from the density of microscopic 
degrees of freedom. We do know that the thermodynamical behaviour of any system 
can be described by an extremum principle for a suitable potential (entropy, free 
energy ...) treated as a functional of appropriate variables (volume, temperature 
,....). If our ideas related to gravitational theories are correct, it must be possible 
to obtain the field equations by extremising a suitably defined entropy functional. 
The fact that null surfaces block information suggests that the entropy functional 
should be closely related to null surfaces in the spacetime. This turns out to be true. 
The Lanczos-Lovelock field equations can indeed be obtained [5] by extremising an 
entropy functional associated with the null vectors in the spacetime. This functional 
depends only on a tensor P"'"='^ which has the symmetries of the curvature tensor and 
is divergence-free in all indices. One can show that the thermodynamic variational 
principle leads to the field equations of a Lanczos-Lovelock model with the same 
entropy tensor. 

In this paper, I will reinterpret these results in a different manner. I will ar- 
gue that it is more natural to project the equations which result from the entropy 
extremisation onto the null surfaces of the spacetime, thereby obtaining what is usu- 
ally called in the literature as the Damour-Navier-Stokes (DNS hereafter) equation 
[SJ[7]- In general this equation is similar to Navier-Stokes (NS) equation but is not 
identical to it because the DNS equation contains a Lie derivative while the stan- 
dard NS equation has a convective derivative. (As I will describe in Section 13.41 
this structural difference exists in all the previous derivations of DNS equation in 
the literature but has not been emphasized). I will argue that it is better to study 
the structure of gravitational field equations in freely falling frames rather than in 
arbitrary coordinates and will show that — in locally inertial frames — the DNS 
equation becomes identical to NS equation. Moreover, the physical interpretation, 
especially aspects related to viscous dissipation, becomes transparent in this frame. 

We thus obtain an entropy extremisation principle to derive DNS equation di- 
rectly, which makes the hydrodynamical analogy with gravity self contained. This is 
in contrast with the conventional approach in which one first obtains the field equa- 
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tions from some (action) principle and then interprets it as a NS like equation. In 
such an approach it seems somewhat mysterious that the equation looks similar to 
NS equation with a thermodynamic structure. 

This paper is organized as follows: I begin in the next section with a short review 
of the derivation of gravitational field equations by extremising the entropy density 
of spacetime. In Section |3.1[ I argue in favour of projecting these equations onto a 
null surface. The technology to describe the extrinsic geometry of the null surface is 
reviewed briefly in in Section 13.21 (A more pedagogical discussion is included in the 
Appendix to make the paper self-contained.) This is used to recast the equations 
in the form of NS equation, with convective derivative, in a boosted inertial frame 
in Section [3731 The DNS equation in a general frame (which has the Lie derivative 
rather than convective derivative) is obtained in Section 13.41 and a comparison of 
these equations and their physical implications are presented in Section 31 Some 
further technical issues related to rescaling of null normal and introducing spacetime 
dependent boosts are discussed in Section 14.11 and the conclusions are presented in 
the last section. This analysis also opens up several further avenues of research which 

1 will briefly mention whenever appropriate. 

We will use the mostly positive signature; the Latin indices, a, 5... go over — 3 in 
the spacetime manifold M while the Greek indices a, (3, ... go over the coordinates in a 
3-dimensional null surface S with signature (—,+,+). When we restrict ourselves to 
the two-dimensional, spatial sub-manifold of this null surface, we will use uppercase 
Latin indices A, B, .... 

2 Entropy density of spacetime and its extremisa- 
tion 

The field equations of Lanczos-Lovelock models were obtained previously [TJ |S] from 
an entropy extremisation principle along the following lines, which we shall rapidly 
review. (The physical motivation and algebraic details are described in detail in 
previous papers, e.g., ref.[I], and will not be repeated here.) The procedure involves 
associating with every null vector in the spacetime an entropy functional Sgrav and 
demanding 5[Sgrav + S matter] = for all null vectors in the spacetime where Smatter 
is the relevant matter entropy. It can be argued, based on the energy flux across 
local Rindler horizons and the associated entropy, that the form of Smatter relevant 
for this purpose can be taken to be 

^matt- / d'^xV^TabU^n" (1) 

Jv 

where n° is a null vector field and Tab is the matter energy-momentum tensor living 
in a general D{> 4) dimensional spacetime. (The integration could be over this 
spacetime volume (k = D) or even over a nontrivial sub-manifold (k < D), say, a 
set of null surfaces. This does not affect the variational principle or the resulting 
equations.) The simplest choice for Sgrav is a quadratic expression in the derivatives 
of the null vector: 

Sgrav ^-4. f (fx^Pab'''Ven''Vdn'' (2) 
JV 

where P^^ is a tensor having the symmetries of curvature tensor and is divergence- 
free in all its indices. (This makes the notation P^^ = Pab^'^ unambiguous and we will 
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often use this placement of indices.) It can be shown that the most general tensor 
which satisfies these criteria is the entropy tensor of a Lanczos-Lovelock theory; 
that is, the P"'""' in Eq. © can be expressed as P"-'""^ = dL/dRabcd where L is 
the Lanczos-Lovelock Lagrangian and Rabcd is the curvature tensor [T] . This choice 
will also ensure that the equations resulting from the entropy extremisation do not 
contain any derivative of the metric which is of higher order than second. (More 
general possibilities exist which we will not discuss in this paper.) The expression 
for the total entropy now becomes: 

5K] = - / d'^x^ (AP^y^n'^Wdn'' ~ Tabu'^n'') , (3) 
Jv 

We now vary the vector field n° in Eq. ^ after adding a Lagrange multiplier function 
A(a:) for imposing the condition Uadn"- — 0. Straight forward algebra (see e.g., 
Section 7.1 of ref. [1]) now shows that the condition S[Sgrav + Smatter] — leads to 
the equations 

i2n^,-T-,+XS^)na^0, (4) 

where TZ'^ = Pb^''^ R'^ijk is the generalization of Ricci tensor to Lanczos-Lovelock 
models. We demand that this extremum condition should hold for all null vector 
fields n". 

At this stage the usual procedure [1, 5' is to use the generalized Bianchi identity 
and the condition V aT§ = to obtain from Eq. Q the equations 

Qt^K-\5tL = h:i + K5t (5) 

where A is a constant and P"'"^'^ — dL/dRabcd- These are precisely the field equations 
for gravity in a theory with Lanczos-Lovelock Lagrangian L. When L (x R, leading 
to P^^ cx ((5^(5^ - S^S'^J, we have oc i?^'^6 « ^'^ ^nd one recovers Einstein's 
equations. 

There exists, however an alternative route which we will explore in this paper, 
that turns out to be more in tune with the entropy extremisation principle and 
emergent perspective of gravity. To do this we note that, while Eq. ^ holds for 
any vector field, once the normalization condition is imposed through the Lagrange 
multiplier, the entropy was originally attributed to null vectors and hence it is 
natural to study Eq. (U) when is the null normal of a null surface S in the 
spacetime and project Eq. (j?]) onto the null surface. This is what we will do. 

In the case of spacelike or tiniclike surfaces with a normal, such a projection is 
straightforward but the null surfaces require a somewhat more careful treatment. 
This is described in the Appendix for those readers who may not be completely fa- 
miliar with the issues involved. In the main text of the paper, we shall use the results 
from this Appendix to proceed further. We will also work with Einstein's theory in 
D ~ 4 for clarity and definiteness and will mention more general possibilities in the 
end. 
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3 Navier- Stokes dynamics of null surfaces 



3.1 Entropy extremum condition projected onto a null sur- 
face 

In a 4-dimensional spacetime manifold Ai with a metric gat, a null surface with 
normal £ will be a 3-dimensional sub-manifold S such that the restriction 7^1, of the 
spacetime metric 1706 to the S is degenerate. It can be easily shown (see Eq. of 
Appendix) that the normal to S satisfies the null geodesic equation £"Va^m = aim 
with a non-affine parametrization, indicated by a non-zero k which we will call the 
surface gravity. Since the normal € to a null surface is also a tangent, this result 
also shows that a null surface can be thought of as "filled by" (a congruence of) 
null geodesies. Using this, we can introduce a natural coordinate system adapted to 
a family of null surfaces in the spacetime. We choose one of the coordinates such 
that = constant correspond to a set of null surfaces with, say, x'^ = on S. 
Let the intersection of S with a = constant surface (Sj) of the spacetime be a 
2-dimensional surface St with transverse coordinates = {x^^x"^) and coordinate 
basis vectors ea = Oa- We will choose the coordinate system such that 

£ = do + v'^dA = da + v^eA; r = (l,t;^,0) (6) 

with with £ • £ = = £ • ea- (That is, we embed 5 in a one-parameter congruence of 
null hypersurfaces corresponding to = constant and choose the other coordinates 
x°,x^ and x^ in such a way that Eq. 1^ holds.) Clearly, £ has the structure of a 
convective derivative if we think of as a transverse velocity field. The line interval, 
which has nine independent functions, has the form in Eq. (I5ip of Appendix. In 
particular the line element on x^ = constant surface is given by; 

ds^ = QAsidx^ - v^dt){dx^ - v^dt) (7) 

The metric on St corresponding to t — constant, x^ = constant is qAB with a well 
defined inverse q^^ . The raising and lowering of the uppercase indices A^B etc. in 
this transverse two-dimensional surface are done using these metrics. 

A more formal way of introducing this metric is as follows; Given a null surface 
with the normal 1°" , we first introduce another null vector fc" with £ • k = —1. (For 
example, in fiat spacetime ii £ = et + is the outgoing null vector, then k could 
be A; = {l/2){et — e^) which is proportional to the ingoing null vector.) The metric 
on the two dimensional surface St orthogonal to this pair is given by the standard 
relations: 

Qab = gab + Lh + hka, Qab^'' = = Qabk'' (8) 

The mixed tensor allows us to project quantities onto St. 

We are interested in the projection of Eq. onto S with n"' — £"■. In contrast 
to timelike or spacelike surfaces, this will now contribute, (i) a term along £ itself 
when we contract Eq. (jj]) with Z° (because £ is also tangent to 5!) as well as (ii) a 
projection to the 2-surface St obtained by contracting Eq. ^ with q\. In both, the 
term involving A will not contribute because £^ = and I'^qab = 0. More formally 
the projection will give (see Eq. of the Appendix) the equations; 

i?mnr^g: = 87rT,„„f'"C; i?„„f "r = 8^r™„r^r (9) 

In other words, the maximization of the entropy associated with the null vectors 
lead to these two equations in a natural, direct manner. 
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Of these, the second equation in Eq. © involvingfJmn^'"^" wih give the famihar 
Raychaudhuri equation (which we will not discuss in detail) while the first one with 
RmnP^fla '^i^ l^ad to an equation that looks like the Navier-Stokes equation — 
which will be of our interest in this paper. 

To obtain the latter in a physically transparent way, we only have to consider 
the A th component (in the transverse direction) of the vector Rmji™' ■ This can be 
expressed in terms of the covariant derivative of ^ by the standard identity 

RraA^ = R^Ai^ = W ^{W - ^^(V^r) (10) 

where the first equality arises from Eq. (Recall that in our notation, Greek letters 
cover 0, 1, 2 while uppercase Latin letters cover 1, 2 in the transverse direction.) To 
rewrite this as a Navier-Stokes equation, we need to re-express the derivatives V 
and V in terms of quantities related to the extrinsic geometry of S. We will now 
briefly overview the concepts involved in the extrinsic geometry of a null surface 
delegating details to the Appendix. 

3.2 Extrinsic geometry of a null surface 

The extrinsic geometry of the null surface is determined by the derivative of the 
normal Vp€ along the tangential directions of S. Because I • V = {l/2)d^t^ = 
the covariant derivative of I along vectors tangent to S is orthogonal to I and hence is 
tangent to S. Therefore V af- is a vector which can be expanded using the coordinate 
basis = on S. (This basis is made of = (9o, Ba) = {(■ — v^eA^eB) = e^.) 
Writing this expansion with a set of coefficients (called Weingarten coefficients) x"' p 
we have 

v„€^xV^-xV^; yj^ = x^^ (11) 

From rVa^'^ = 1°"^ = J"" = ^^^^ , it follows that ^" is an eigenvector of 
the matrix with the eigenvalue k. The coefficients determine the extrinsic 
geometry of the null surface we are interested in. We will now study the different 
components of the 3x3 matrix x^a- 

To do this, we define the quantity = ^^^^x^ a components of which are 
given by 

= 7M/3X"a = ea ' ^ 0^ = V^(£ • e„) - £ • V^e„ = -l.V'^f, (12) 

where we have used £ • = 0. This quantity is clearly symmetric (0q^ = 0^q) 
because Christoffel symbols are symmetric. Further, from the result Ih^ — 0, 
we have the identity Qafi^^ — which, on expansion, using If' = (1 , TJ"^), gives 
0qO = —QaBV^. This implies Qqa — ©ao = —Qabv^ which, in turn, leads to 
000 = —'doBV^ = v^v^Qab- Thus all components of Qai3 can be determined in 
terms of the three components of 8 as and the metric coefficients v^. 

The components &ap have a direct geometrical significance as the Lie deriva- 
tive along £ of the metric qij on S. Using the restriction of the 4-dimensional Lie 
derivative of gij to 5t, we get 

Qab = ^££qAB (13) 

in the adapted coordinate system. This structure is similar to the standard formula 
K^i, = —{l/2)£jih^u for the extrinsic curvature oi at = constant spacelike surface 
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in the (1 + 3) decomposition. In fact, evaluating Eq. ((T3|) in the coordinate system 
adapted to S we get a result very similar to the familiar one for K^^, (see e.g., Eq. 
(12.21) of ref. [5]): 



1 ^ dqAB \ 



Qab = 2 y^AVB + Dbva + (14) 

When qAB is independent of t, we see that Qab is essentially the shear tensor of a 
velocity field va ■ This is the key reason why the projected equations can be interpreted 
as the Navier-Stokes equation. 

We will also need to define three more quantities uja = (i^OjI^a) in terms of the 
Weingarten coefficients by 

= x\ = v„£° = r^"„£^" (15) 

The utility of Qab and {ujo,uja) arises from the fact that we can determine all the 
components of x°' p from QaBt^u and the metric coefficients. Among the compo- 
nents of p already have the direct relations x^o = '^OiX^a — ^A- As for the 
remaining components x'^o-X^b- note that: 6ob = X^alt^B = x"o7ob + X^oIab 
giving 

X\ = (eoB - x"o70B)g^^ - -{Qbcv'' - LOoVB)q''^ (16) 

Similarly we have 

X^A = - x\lOB)q''^ = i^AB + iOAVB)q''^ (17) 

So the covariant derivative of the normal vector t — and thus the extrinsic 
geometry of 5 — is completely characterized by the set {Qab — —^nJ'^AB^^o ~ 
P^T^Q, UJA — ^'"rmyi)- Confined to the transverse components, Eq. ([TT]) also gives 
the relation: 

Va^b = XAB = [6 AS + wbwa] (18) 

which will turn out to be useful. 

Being a symmetric tensor, @ab can be expressed in terms of its irreducible parts, 
viz. the trace 9 = 0^ and the trace-free part cr^ by 

+ {l/2)5ie (19) 

In the adapted coordinates, 6 is given by the expression 

e^^^\n^+DAV^ (20) 
and is related to the divergence of i through the relation 

Vat = V^r = Va^^ + Vo^° = {0 + VALo^) +LOo = e + = 6 + k (21) 

In arriving at the third equality we have used Eq. (jl8p and in arriving at the first 
and the fourth equality we have used Eq. The last relation arises because 

X^/f" = k£° = K^ UJ^^^'. 
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3.3 The Navier-Stokes equation in the boosted inertial frame 

Using these results, we are in a position to rewrite Eq. ^TU\\ in terms of the variables 
related to the extrinsic geometry of the null surface. Using Eq. (fTTI) and Eq. (PTjl we 
can rewrite Eq. (fTO)) as: 

" = R^A^^ ^W^x%-dA{i^ + 0) (22) 

Let us consider this expression around any given event in a frame in which the 
Christoffel symbols vanish but their derivatives do not. We will assume that the 
metric coefficients become constant and have their usual diagonal form except for 
keeping ^ but a constant. (That is, we are working with a Lorentz-boosted 
local inertial frame; we will comment on this choice later on in Sec. 14.11 ) In that 
case, we can write V^x^'a in Eq. (j22l) as: 

V^X% = 5^X% = 5oX% + 9sx''a = 9oc^a + 5b(<z^^Xca) 

= doujA + Ob [q""^ [QcA + i^AVc)] = doiOA + asef + v^'Obuja (23) 

Each equality is obtained by discarding terms involving Christoffel symbols but 
retaining terms containing their derivatives; we have also used Eq. (I15p and Eq. (|18p . 
(Note, for example, that is proportional to Christoffel symbols and hence vanishes 
but its derivative need not.) Thus we get 

RmAi"" = (do + v^'ObHa + de (<J^ + ^eS^] - dA{K + 9) 



= {do+v''dB)iOA + dBCJ^-dA[n+h^ (24) 

where we have used Eq. (I19p . The first equation in Eq. ([9]) now becomes: 

(So + v''dB)uJA + ObTa -dA(^fi+ = SttT^^^™ (25) 

Again note that vanishes in the local inertial frame but its derivative ObCa '^i^l 
not. Rewriting Eq. (j25p in the form 

id. + ^'^b) (-^) = - ^^dAe - OA (^) - T^Ai- (26) 

we see that Eq. (|26)) has the exact form of a Navier-Stokes equation for a fluid with (i) 
momentum density —uja/Stt, (ii) pressure {k/&t:), (iii) shear viscosity coefficient rj — 
(I/IGtt). (Note that in the conventional NS equation the viscous tensor 2riag+S,6^6 
is defined with an extra factor 2 for shear viscosity.) (iv) bulk viscosity coefficient 
^ = — l/167r and (v) an external force Fa = T^a^"''- The NS will also have term 9uja 
which vanishes in the local frame because 9 = but, of course, we can formally add 
this term (which is numerically zero) to Eq. (j26p to complete the structure; but the 
derivative OaO is non-zero, which allows us to uniquely determine the bulk viscosity 
term. 

For this interpretation to be strictly valid, it is also necessary that cr^ has the 
form of a trace-free shear tensor built from the velocity field va- From Eq. (jl4p we 
see that this is true provided [dqAB/dt] = when we have: 

ds'J^ = Idsiid^VA + d^B) - (ac«^)] (27) 
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In the literature, Eq. (|26l) or its analogues are often called as Navier-Stokes equation 
even when (dqAB / dt) is non-zero. This is — strictly speaking — incorrect especially 
because the metric qab has no fluid dynamical interpretation in general. We shall 
come back this point in the next section. 

It should also be noted that, in the conventional NS equation, the momentum 
density that is transported by the fluid is usually coUinear with velocity. In our 
case, the uja and va are in general unrelated and — in fact — the physical meaning 
of LUA is unclear for a general null surface. (It does have a physical meaning in 
the context of the event horizon of Kerr black hole, for which the formalism was 
originally delevoped; in that context, uja can be related to the angular momentum 
of the black hole.) 

It is clear that 77/^ — —1 also makes any interpretation in terms of standard 
viscous dissipation unnatural for a general null surface In the literature, this 
interpretation has been provided in situations when teleological evolution is accept- 
able especially in the context of black hole event horizon. More importantly, the 
ratio between shear viscosity and entropy density s = 1/4 of the null surface is given 
by 77/s = l/47r, when we use the fact that any null surface can be attributed an area 
density of entropy which is (1/4) in Einstein's theory. Similar results are known in 
much more complicated situations arising from string theory but they presumably 
have the basis in Eq. 

For completeness we will also state the Raychaudhuri equation arising from the 
second equation in Eq. ^ which, in the boosted local inertial frame reads (see 
Eq. (p5|) of Appendix and note that 6„m = = 6* in this frame): 

Rmni'^r = -{do + v^'ObW = 87rr™„^"r. (28) 

which relates the evolution of 6 to the energy flux across the null surface. In the 
context of local Rindler frame, the right hand side of Eq. can be interpreted 
as proportional to the heat flux making this an entropy balance equation. We shall 
briefly comment on this interpretation in Sec. 15.11 

3.4 The Damour-Navier-Stokes equation in arbitrary frame 

The entire analysis analysis can be repeated in an arbitrary coordinate system to 
obtain a generalized Navier-Stokes equation, first obtained by Damour [6]. This is 
derived in the Appendix (see Eq. (jTO]) ) and has the form: 

where the various terms are defined below. The metric is given by Eq. ([5]); the 
rim is the projection ftm = Qm^n = — nkn where a;„ = PVjkn] note that, 
since kA = 0, we have i^A = as far as the transverse components, which we are 
interested in, are concerned. The Qmn is the projection of the covariant derivative: 
©mn = QmQn^a^b- The Di is the covariaut derivative operator defined using the 
metric qab in the two-dimensional subspace, while £(, is the Lie derivative with 
respect to t. 

A comparison of Eq. with Eq. shows that — while most terms have 
an one-to-one correspondence, with suitable projections replacing the operation of 
taking transverse components — there are two crucial differences which we will now 
discuss: 
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(i) In the literature one usually sees the interpretation of as the viscous 
shear tensor of a fluid with a velocity field Va without bothering about its internal 
structure. But, as we have already mentioned, this tensor contains an extra term 
involving the time derivative of the transverse metric, (dqab/dt) which has no direct 
fluid dynamical interpretation. Similarly, the bulk viscosity term involving 9 contains 
the time derivative (d/dt) In y/q. (See Eq. pi)) and Eq. (1^ : of course, (dqab/dt) = 
in the frame we are working with, but its derivative dmidqab/dt) does not vanish, 
which is what leads to the difficulty) . To get out of this difficulty, one needs to make 
certain choices for the metric and coordinate system so that Qab is independent of 
t. (In some other contexts one can avoid this term by rescaling £ or by working in 
a perturbative series for the metric; we shall briefly comment on these possibilities 
in Sec. 14.11 ) There is very little discussion of this issue in the literature and the 
term involving a™ is usually called shear viscosity term without worrying about its 
internal structure. 

(ii) The second, and probably more crucial, difference is the following: The DNS 
equation in a general coordinate system, given by Eq. ((29)) . contains the Lie derivative 
with respect to £ while the standard fluid dynamical, non-relativistic, NS equation 
contains the convective derivative with respect to £. This difference can be explicitly 
verified by expanding out Eq. ()29l) in the coordinate system adapted to the null 
surface in which £ ~ + v^Da (see Eq. (|75)) in Appendix): 

{do + v^Db)uja + LOB Dav^ + Oloa + Db<Ja -Da(^k + ^9^ = SnTmni'^qA (30) 

The key difference is second term in the left hand side, ujb Dav^ , involving the 
derivative of the velocity field which is absent in the standard NS equation, making 
the two equations structurally different. (This term vanishes in the local inertial frame 
since uja vanishes in that frame, making the structure of DNS equation identical to 
NS equation.) More formally, the convective derivative term Dtfta = q^P^^ m^n 
defined with correct projections, is related to Lie derivative term by: 

Dt^la = Cr^V„r!„ = qZ'£e^n - Qa^m (31) 

Using this we can rewrite Eq. (|^^ as: 

Dt^a + e^r!,„ + 9na-Da(^K+'^^+ D„-,a'^ = 87rT,„„ (32) 

which clearly brings in an extra term, which is the second term in the left hand side. 

Our Eq. ()29p and Eq. (|5(I)) agree with the original results obtained by Damour 
[see e.g., Eq. (1.55) of ref. [5] or Eq. (15b) of ref [7]]. In ref. [B], Damour also 
quotes the standard NS equation (see Eq. (1.47) of ref [7] ) which, of course, differs 
by the term involving the derivative of the velocity from the Eq. (1.55) of the same 
work. Similar difference exists in ref. [10] between the DNS equation [Eq.2.67] and 
the NS equation [Eq. 2.68]. This result also has been obtained in the membrane 
paradigm approach in ref. [TT] and our Eq. matches with Eq. (2.14) of ref. [TT] : 
the second term in the left hand side of their Eq (2.14) is the extra term Q^flm 
of Eq. making it different from the NS equation. (Price and Thorne, however, 
does not call this equation Navier-Stokes equation, probably because of the extra 
term, and refers to it as Hajicck equation.) Since our derivation in the Appendix is 
modeled exactly as in ref. [12j . it, of course, matches with their equations in Section 
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(6.3). This review, as well as many other papers in the literature (like, for e.g., 
call our Eq. ([23) or Eq. ([50)) the DNS equation even though there is an extra term. 
It is true that the Lie transport arises more naturally than convective transport in 
curved spacetime and one may think of it as a natural generalisation; but the extra 
term in a PDE makes its structure (and solutions) very different and hence a careful 
distinction needs to be made. 

We have, however, shown that one can indeed think of DNS equation as identical 
to the NS equation provided we view it in the (boosted) local inertial frame. Since 
u)a is proportional to the Christoffel symbols, the extra term vanishes in the local 
inertial frame making this interpretation easy. This procedure actually has a broader 
domain of relevance which I will now describe. 

4 Fluid dynamics of null surface and the issue of 
viscous dissipation 

In any generally covariant theory, one tries to express the equations in a form valid 
in arbitrary coordinate systems in order to exhibit the symmetries explicitly. While 
this is certainly important to check the diffeomorphism invariance of the theory, it 
has two shortcomings in certain contexts: 

(a) The metric coefficients will depend on the spacetime coordinates both in 
curvilinear coordinates as well as in a curved spacetime. Principle of equivalence 
implies that one cannot meaningfully ask "how much" of this dependence arises 
due to the choice of the coordinates and "how much" arises due to the curvature 
of spacetime. When we are interested in the structural aspects of the theory, the 
spacetime dependence of the metric due to the choice of coordinates acts as an extra 
baggage which one would like to avoid. Working in the local inertial frame allows 
us to do this. 

(b) The emergent paradigm of gravity suggests that one needs to accept an 
intrinsic observer dependence in all thermodynamical variables. In particular, one 
can attribute an entropy, temperature and related dynamics to local Rindler horizon 
as perceived by the corresponding Rindler observer in the spacetime even though a 
freely falling observer will not attribute any thermodynamical features to the same 
null surface ([HJ [TS]; also see Sec. 4.4 of ref. [T]). But, while we use the observer 
dependent thermodynamics to obtain the field equations of gravity in this approach, 
the final result should be (and is) independent of the Rindler observers and must be 
generally covariant. Therefore while projecting the field equations to a null surface 
and interpreting them as NS equation of a fluid living on that null surface, it is 
crucial to ask which part of the field equation survives in the freely falling frame. 
Our analysis shows that the terms involving (i) the derivatives of the stress tensor, 
(ii) the derivative of the pressure and the (iii) external momentum flux term remain 
non-zero in the freely falling frame because they involve derivatives of the Christoffel 
symbols. This is in spite of the fact that the viscous stress tensor and the pressure 
themselves vanish in the freely falling frame. What is relevant for the equation to be 
interpreted as the NS equation is the existence of a viscous force on the fluid, arising 
from the gradient of the viscous tensor (and pressure) , rather than the viscous tensor 
(or pressure) itself. Of course, this is very counterintuitive compared to normal fluid 
mechanics. If flow of water exhibits viscosity, it will also have non-zero viscous tensor 
in the same frame of reference! But here, we are merely calling some combination of 
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Christoffel symbols as viscosity tensor (see Eq. and their derivatives involving 

certain combinations of Ricci tensor as viscous force. Obviously, the former can 
vanish without the latter vanishing in a local incrtial frame. 

This feature has important implications for the characterization of viscous dis- 
sipation in the current context. In the conventional fluid dynamics based on NS 
equation, the viscous dissipation will be proportional to terms involving aAB^^^^ 
and 6^ and — in usual fluid mechanics — (t^bct"*^, and dscr^^ will all be nonzero. 
In our case uab and 6 vanish in the freely falling frame and the inertial observers 
in spacetime will not see any dissipation. (The Raychaudhuri equation in arbitrary 
coordinates — Eq. ((65)) in the Appendix — has these quadratic term in the viscous 
tensor which disappear in Eq. (|28|) valid in the locally inertial frame.) This is reassur- 
ing since we do not probably expect a continuous, observer independent, dissipation 
to take place in spacetime. But the force on the viscous fluid, which depends on the 
gradient ObO'a and OaO, does not vanish even in the freely falling frame showing that 
this force has an observer independent existence. It might appear, at first, paradox- 
ical that a viscous force term exists for the fluid but no viscous dissipation. But, as 
we said above, algebraically this is no more paradoxical than the fact that effects due 
to curvature involving derivatives of Christoffel symbols can be present locally even 
when the Christoffel symbols vanish at a point. The result obtained above is just a 
translation of this well-known fact in the language of null surface dynamics. In fact 
the real paradox would be if there is an observer independent viscous dissipation in 
spacetime! 

There have been attempts in the literature to interpret the equations for the 
fluid on the null surface, especially the Raychaudhuri equation, as describing entropy 
production ([HIIZIITT]; for more recent work based on emergent approach to gravity, 
see e.g., [16]). The early work (e.g. [6l[7j[TT] leads to an entropy production that is 
acausal (or teleological, as it is often called) involving the entire future history of the 
null surface. In this context, we need to distinguish between two different physical 
situations. This teleological feature is probably acceptable in the context in which 
DNS equation was originally derived and applied, viz., to black hole event horizons. 
Event horizons are in any case teleological even in their definition and hence this is 
not a surprise. 

But the DNS equation is applicable to any null surface in the spacetime, including 
patches of the local Rindler horizons. In the emergent paradigm of gravity one relies 
heavily on local Rindler horizons which are ordinary null surfaces in spacetime and 
not event horizons corresponding to some specific solutions of field equations. In fact, 
to obtain the field equations from a suitably defined entropy density of spacetime, one 
uses the local Rindler horizons as off-shell constructs in the theory and hence they 
cannot be linked to horizons arising in on-shell solutions. But since, in our derivation 
of DNS we never assigned any special property to the null surface in question, the 
resulting equations are applicable as equally to a local Rindler horizon in a given 
spacetime as to an event horizon of a black hole obtained as specific solution. This is 
why we need to carefully distinguish generic features of interpretation which are valid 
for any null surface from some of the features which might have limited validity in 
the specific context of black hole horizon. The dissipation and generation of entropy 
with accompanying irreversible thermodynamics may be acceptable in the case of 
physical processes involving black hole horizons, say, but such dissipation is difficult 
to interpret in terms of local Rindler horizons. So, again, it is probably welcome that 
o'abO'^^ and 6^ vanish in the freely- falling frame but the NS equation for the fluid 
remains valid. From this point of view, it seems better to interpret the equations in 



12 



the freely falling frame. 

Incidentally, it does not seem to have been emphasized in the literature that 
Einstein's equations have a very simple form in the freely falling frame. Using the 
exact expression for Ricci scalar 

Rab = -^drniV^ T™,,) - d^^V^) " r'"„,r" (33) 

we find that, in the locally inertial frame Einstein's equations take the form: 

Rab = dr,X'\b - dadbV^ = ^Tab " \riabT) (34) 

In fact projecting Eq. p4p to a null surface (when riabT term will not contribute) we 
can re-derive our NS equation. Since all the terms in the left hand side of Eq. 
correspond to specific combinations of curvature tensor components we see that the 
DNS equation in the boosted inertial frame essentially provides mapping between 
curvature components and the relevant fluid variables. This aspect will be explored 
in a separate publication |17) . 



4.1 Spacetime dependent boosts and rescaling of the null nor- 
mal 

In obtaining Eq. (j25p we have purposely used a boosted local inertial frame in which 
all the metric coefficients take the standard Cartesian form except goA = —va- One 
could have, of course, repeated the derivation with va [and its derivative which, 
being proportional to Christoffel symbols, are anyway zero] set to zero but with 
the second derivatives of va (which contribute to the viscous force term in 9b cr^) 
remaining non-zero. We kept the va non-zero in the original derivation just to 
exhibit the nature of convective derivative and to contrast it with the Lie derivative. 
Physically, one can think of non-zero va as arising because of a transverse velocity 
for the observers (e.g., local Rindler observers with a constant velocity drift parallel 
to the Rindler horizon). 

The derivation of DNS equation does not pre-suppose any special form of the 
metric. In certain contexts, one can choose the metric as well as the normalization 
of (. in such a way that the extra term in the Lie derivative is cancelled out. This 
arises because the different terms in the DNS equation scale differently when we 
rescale the null vector thy & spacetime dependent factor 7(0::). (Being a null vector, 
£ has no natural normalization, without us making some additional assumptions.) 
This transformation changes the components of t to the form (7, 0, 7"^^) suggesting 
the interpretation of a spacetime dependent Lorentz boost with a 7 factor. One can 
easily show that when t I' ~ "ft the Qab^ ^ab and 9 scales by a factor 7 and the 
vector k goes to "f~^k. But n and flA change inhomogeneously as: 

k' = 7(k + ln7); fl'j^ — Ha + DaIhj (35) 

Because of Eq. ([55)1 . the DNS equation does not retain its form under such a rescaling. 
One can, by hand, choose the value of 7 to cancel out the extra term in the Lie 
derivative thereby reducing the equations to the standard N-S form. (This has been 
done, for example, in ref. [13] with a different motivation.) 

A similar mathematical situation arises in a more general context of a set of 
transformations which are usually used in the literature. Consider a metric written 
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in a coordinate system with goa = so that the line element is 

ds^ = -F^dt^ + U/sdx^dx^ (36) 

We will now substitute in this metric dt = Uadx"" where Ua has the components 
M° = 7(1, u"). If we take to be constant with 7 = (1 — u'^)^^/^, then such a trans- 
formation is equivalent to a Lorentz boost and, of course, leaves the physics invari- 
ant. But when and 7 are functions of spacetime coordinates, the transformation 
changes the character of geometry. (Without additional constraints, dt = Ua{x)dx°' 
will not even be integrable and dt will not be an exact differential; however, one can 
certainly use this "rule" to modify the metric.) Straightforward algebra now shows 
that the resulting metric can be expressed in the standard (1+3) form as: 

ds^ = -N^dt^ + Kfs{dx°' + N"dt){dx'^ + Nl^dt) (37) 

where 

^2 ^ fV(1-J^'7V)"'; h^p ^ Up - F^-f%^vp 

= FV^'^(1 -^^7^^'^)"^ N^^F^-f^Vf, (38) 

Obviously, such a modification of the metric is not a symmetry of the theory when 
w°(x) are spacetime dependent functions; for example, if the original metric in 
Eq. pop satisfied Einstein's equations with some source, the modified metric in 
Eq. ((57)1 will not be a solution unless very special conditions are satisfied. We will 
provide a general discussion of these transformations in a future publication [17] but 
will point out one simple feature here which is relevant to our discussion. 

Notice that if we consider the lowest order change, linear in the boost velocity, it 
is essentially the generation of the cross term with N^^, « F'^v^. This has the effect 
of generating an extra term to the extrinsic curvature of the t = constant surface 
given by 

6K^p = ^[Do.iF\p) + Dp{F''v^)] (39) 

where D is the covariant derivative based on the spatial 3-metric. When the original 
coordinate system was synchronous with F = 1, this extra term has the form 

6Kap = ^ [D^vp + Dpva] (40) 

Further, if we treat one of the spatial coordinates, say, x^ differently from the two 
transverse coordinates x^ = {x^,x^), then the Qab of the null surfaces orthogonal 
to null vectors in x'^ — x^ plane will pick up a term of the form SQab = (1/2)[-Da'1's + 
Dbva] (see Eq. (1611) in the Appendix). This, in turn, corresponds to the addition of 
an extra term 

S<TAB = ^[{Davb + Dbva) " qAB{Dcv^)] (41) 

to the trace-free part of Qab which has the form of the shear viscosity tensor. This 
analysis shows how a drift velocity in the transverse direction to a null surface can 
lead to the viscous stress tensor of the correct form. Several discussions of this topic 
in the literature are essentially special cases of the transformation in Eq. ([55)) with 
very specific choices for F, fap and Ua- One is led to an NS equation in all these 
contexts essentially because of the result in Eq. (|1T|) . 
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As a simple illustration of the same effect, consider a metric of the form 

ds^ = -N^dt^ + T7T + SABidx"^ - v^dt){dx^ - v^dt) (42) 

which could represent the local Rindler frame when = Q and N{x) cx x with 
the surface a; = being null and acting as the local Rindler horizon. Since such 
a metric represents flat spacetime, it obviously satisfies source-free Einstein's equa- 
tion. If Tj"^ 7^ hut constant, the metric represents a boosted frame in which the 
Rindler observer has a constant drift velocity in the transverse direction. When 
_ yA^^^^A-^ ^jjg metric is still static but will not, in general, satisfy source-free 
Einstein's equations. This introduces a spacetime dependent boost velocity with 
non-zero shear. For a general N[x) the null surfaces with translational symmetry in 
the y, z directions are given by /(t, x) — constant where 

( dx , , 

/(i,^) = -i + y^ (43) 

The normal to this null surface can be taken to be ia — iV^(— 1, A^^^, 0, 0). (The 
overall normalisation, as we have discussed, is not unique but this form is based on 
Eq. ([5^ of the Appendix, which has some level of naturalness.) A direct calculation 
shows that, for these null surfaces, we have 

Sab = \^dAVB - Qb^a^ (44) 

with il^ = {Xl'^)dxVA and k = dxN^ ■ It is obvious that Qab has the form of a shear 
tensor for a fluid with velocity v'^ . So the null surface fluid equations will pick up a 
viscous tensor arising from this shear. This is essentially what happens even in more 
general contexts. 

Finally, it may be noted that the validity of DNS equation on a null surface forms 
the basis for the structural similarity noticed between gravity and fluid mechanics 
in the context of string theory motivated approaches in AdS etc (see e.g., [H]). 
Once the field equation is reduced to the DNS form (or its generalization in a higher 
dimensional AdS context), one can attempt to solve the equation — which is an 
issue we have not addressed — and obtain the metric systematically using a long 
wavelength approximation of the fluid mechanics. This would involve expanding the 
metric in a series gab = gab° +^gab^ +£^gab^ ■■■■ indexed by a book-keeping parameter e 
and, simultaneously introducing a derivative expansion by rescaling the coordinates 
by X* — > ex^ in the argument. This will allow one to solve for the metric order by 
order in the long wavelength approximation exactly as in the case of fluid mechanics. 

5 Conclusions 

The similarity between gravitational fleld equations and the Navier-Stokes equation 
in the context of null surfaces was known (at least to the relativity community) for 
several decades now. However, it had not found any specific utility except possibly as 
a conceptual tool in membrane paradigm until recently when similar results, inspired 
by string theory, attracted attention to this subject. In this backdrop, the results of 
this paper have the following significance and implications: 
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5.1 Structure of spacetime entropy density functional 

In all the previous approaches to the DNS equation, one starts from the gravitational 
field equations (which would have been obtained by some standard field theoretic 
procedure), projects them onto a null surface and reinterprets the resulting equation 
as the DNS equation. This is somewhat unsatisfactory conceptually since one would 
have expected a fluid dynamical equation to arise from an entropy maximization 
principle or something similar. Given the fact that Einstein's field equation itself 
has clear thermodynamical interpretation, it is obvious that DNS equation itself must 
be derivable in a direct manner from entropic considerations. We have achieved this 
goal in this paper by reinterpreting the previous result [Tl [5] of obtaining the field 
equations from extremising entropy density of spacetime. 

It is conceptually satisfying that the extremum condition for the entropy density 
of spacetime (directly related to null surfaces because they can act as local Rindler 
horizons) given by Eq. Q has a natural interpretation when projected to a null 
surface. As Eq. ([9]) shows, such a projection — involving contractions with either 
or — makes the Lagrange multiplier term in Eq. 1^ disappear and gives a direct 
relation between Ricci tensor and energy momentum tensor. 

It may also be noted that the gravitational entropy density — which is the 
integrand Sgrav oc {— P^^V ci"'^ d^) m Eq. (0) — obeys the relation: 

(X {-P^J^d') « (V„r - 5y,P) (45) 

where the second relation is for Einstein's theory. This term is analogous to the more 
familiar object = K^ — 6^K (where Kab is the extrinsic curvature) that arises in the 
(1+3) separation of Einstein's equations. (More precisely, the appropriate projection 
to 3-space leads to f^.) This combination can be interpreted as a surface energy 
momentum tensor in the context of membrane paradigm because tab couples to Sh"^^ 
on the boundary surface when we vary the gravitational action ( see, e.g., eq.(12.109) 
of [5]). In fact, one obtains the results for null surfaces as a limiting process from the 
time-like surfaces (usually called stretched horizon) in the of membrane paradigm 
[TT| . Equation (H5|) shows that the entropy density of spacetime is directly related 
to and its counterpart in the case of null surface. 

This term also has the interpretation as the canonical momentum conjugate to 
the spatial metric in (1+3) context and Eq. P5)) shows that the entropy density leads 
to a similar structure. This will form the basis for generalising the DNS equation to 
Lanczos-Lovelock models in a future publication '17| . 

Further, the functional derivative of the gravitational entropy in Eq. ^ has the 
form, in any Lanczos-Lovelock model: 

^-^ocTZabe'ocJa (46) 

Previous work [1] [3l [M] has shown that the current Ja = 2TZab^'' plays a crucial 
role in interpreting gravitational field equations as entropy balanced equations. In 
the context of local Rindler frames, when arises as a limit of the time-like Killing 
vector in the local Rindler frame, Ja can be interpreted as the Noether (entropy) 
current associated with the null surface. In that case, the generalization of the two 
projected equations in Eq. ([9]) to Lanczos-Lovelock model will read as 

Jar^^Tabte"; Jaq^ra = \Tabtql (47) 
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which relate the gravitational entropy density and flux to matter energy density and 
momentum flux. The second equation in the above set becomes the DNS equation 
in the context of Einstein's theory. All these results, including the DNS equation, 
will have direct generalization to Lanczos-Lovelock models which can be structured 
using the above concepts pTJ . We again see that all these ideas find a natural home 
in the emergent paradigm. 

5.2 Comparison of DNS and NS equations 

We have used this occasion to clarify some issues related to the interpretation of 
DNS equation as describing a viscous fluid. As noted in Section 13.41 there are, in 
general, two crucial differences between DNS equation and standard NS equation. 
The first is the appearance of time derivatives of the transverse metric qab in the 
definition of shear viscosity tensor because of which it cannot be interpreted properly 
in general. In fact this term has no hydrodynamical interpretation except in very 
specific contexts or in a perturbative series. One can, however, avoid this difficulty 
by choosing specific class of metrics and coordinate systems. 

The second difference has to do with the appearance of Lie derivative rather than 
the convective derivative in the momentum transport equation. As I argued, this 
requires one to work in the boosted local inertial frame for proper interpretation. 
Somewhat surprisingly, these two differences have not been explicitly discussed in 
published literature, to the extent I know. 

5.3 Dissipation: a new level of observer dependence 

The existence of viscosity in the fluid equation raises questions related to possible 
dissipational effects. I am not comfortable with the notion of continuous, observer 
independent, dissipation in the spacetime and have stressed the fact that the dissi- 
pational terms involving (Tabcr"'^ and 6^ vanish in the boosted inertial frame which 
I consider to be appropriate for the interpretation of DNS equation, in view of the 
comments in Sec. 13.41 and Sec. 15.21 As pointed out earlier, in such a frame, the 
derivatives of the viscous tensor (related to spacetime curvature) are non-zero while 
the viscous tensor (proportional to Christoffel symbols) vanish. Of course, an ob- 
server who is not freely falling will perceive non-zero dissipational terms while a 
freely falling observer will not. This may sound paradoxical but observer depen- 
dence in the thermodynamic description of horizons is very well known and need not 
cause any (new) surprise. 

A general framework incorporating the observer dependence of thermodynamic 
variables and providing a translation table for physical phenomena perceived by 
observers in different states of motion is currently lacking. This aspect deserves 
further investigation. 
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Appendix: Null surfaces 



Consider a 4-dimensional spacetime manifold A4 with a metric gab- A null surface 
in this spacetime is a 3-dimensional sub-manifold S such that the restriction ^^ly of 
the spacetime metric gab to S is degenerate, i.e., there exist vectors in S such 
that 7^1/ = 0. (Recall that we use the mostly positive signature; the Latin indices, 
a, b... go over — 3 in while the Greek indices a, (3, ... go over the coordinates in 
a 3-dimensional null surface S with signature (— , +). When we restrict ourselves 
to the two-dimensional, spatial sub-manifold of this null surface, we will use upper- 
case Latin indices A, B, .... Our results are easily generalizable to a I?— dimensional 
manifold but we will stick with 4 dimensions for notational simplicity.) 

The normal to 5 is a null vector field in the spacetime which can be written 
in the form £a = fda4> where / and are scalars in the spacetime and cf) is constant 
on S. Straightforward algebra using the fact ViV = VjVicj) shows that 

r Va^™ = {td. In /) C = (48) 

Therefore the rescaled vector L™ = f^^£"^ is null and satisfies the geodesic equation 
with afRne parameterization; i.e., L^V^L^ — 0. Hence a null surface can be thought 
of as "filled by" (a congruence of) null geodesies. Equation (|48)) shows that also 
satisfies the geodesic equation but with a non-affine parameterization if k is non-zero. 
When the null surface corresponds to a black hole horizon in an asymptotically flat 
spacetime, there is a natural choice for such that k can be identified with the 
surface gravity of the black hole horizon. We shall continue to use this terminology 
and call k in Eq. (|48)) as the surface gravity of S even when it depends on the choice 
of normalization for the null vector field. In the case of S being a local Rindler 
horizon, one can relate k to the acceleration of the congruence of Rindler observers 
which one is considering; the arbitrariness in the normalization translates into the 
arbitrariness in the choice of the acceleration for the Rindler observers. 

The fact that the null surface is spanned by the null geodesies allows us to 
introduce a natural coordinate system adapted to a family of null surfaces in the 
spacetime as follows: We choose one of the coordinates such that = constant 
correspond to a set of null surfaces with, say, a;'^ = on S. Let the intersection of 
S with a a;° = constant surface (Et) of the spacetime be a 2-dimensional surface St 
with coordinates = (2:^,2;^) and coordinate basis vectors = 9a- At any point 
V in Et there will be one future pointing null direction orthogonal to the 2-surface 
St. We choose £ ai V to be in this direction with £ • ba — on St- We can now 
erect a coordinate system in the neighborhood of Et by the following choice: (a) 
The coordinates x^,x'^, x^ are taken to be constant along the geodesies starting from 
each point V{x^,x^,x^) in Ej in the direction of £"-{x); (b) x'^ is chosen to be the 
affine parameter distance along these geodesies with A = t on Et . 

In such a coordinate system, £ = d/dx^ so that = 5^ and the condition £^ = 
translates to goo = 0. Further, the geodesic condition with affine parametrization 
gives 

= t^,t - Vo^ = = ]^g'^\2dogab - ^6500) = g^'^dmb (49) 

requiring 9o5o& — along each geodesic. But initially, on Et we have 70A ^ £-eA = 
and 700 = because £ is orthogonal to the basis vectors ba as well as to itself. This 
requires (7ooi7oa) to vanish all along the geodesic. The line element will now take 
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the form: 



ds^ = -N'^dt^ + (j^dx^ + eNdtj + qAB{dx^ + m^dx^){dx^ + dx^) (50) 

with e = ±1. The x^ = constant surfaces are nuU with the hne element ds^ — 
qAsdx^dx^ because in the (degenerate) metric 7^^ (with /i, = 0,1,2) the coef- 
ficients 7o^ vanish: i.e., 700 = 7oA = 0. For example, when N = M = l;qAB = 
Sab ^T-^ = we recover the usual null coordinates x^ = z + t or x^ ^ z — t for the 
two choices of e. 

In this construction of coordinates, we have embedded 5 in a one-parameter 
congruence of null hypersurfaces corresponding to = constant and the other 
coordinates x'^,x^ and are constructed in such a way that {. = d/dx^ is a null 
geodesic field in the neighborhood of S. While this is always possible it is often 
advantageous to use a slightly less restrictive coordinate system in which the geodesic 
condition is relaxed but the coordinate x^ is constant on the null hypersurfaces of 
interest. This can be achieved by allowing for non-zero 700 and ^oa- In such a case, 
the coordinates can be chosen such that the line interval becomes: 

ds^ = -N'^df+(^dx^ + eNdt^ +qAB[dx^ ~v^dt+'m^dx^){dx^ -v^ dt+m^ dx^) 

(51) 

The metric on St corresponding to < = constant, x^ ~ constant is qAB with a well 
defined inverse q"^^ . The raising and lowering of the uppercase indices A,B etc. in 
this transverse two-dimensional surface are done using these metrics. We also have 
= M(det qAB)^^^- The null vector field will now he £ = da + v^dA ~ da + v^eA 
with components £° = (1, t''^, 0) with f. • £ = = £ • ba- Clearly, £ has the structure 
of a convective derivative if we think of as a transverse velocity field. 

Prom the form of the metric in Eq. (|5ip it can be explicitly verified that ^^u^^ = 
where the Greek letters run over 0, 1,2 on the null surface S with x^ = constant. 
Therefore the 3-metric on S is indeed degenerate. This fact requires us to be careful 
in characterizing the extrinsic geometry of S since raising and lowering of indices will 
not be possible with a degenerate metric. As outlined in the main text, this can be 
done by expanding V a£ using the coordinate basis = on S. For many purposes 
it is convenient to have a more formal approach in terms of tensorial objects defined 
in the 4-dimensional spacetime and their projections onto St ■ We shall now indicate 
how this can be done along the lines described in the excellent review, ref. |12| . 

We will begin by introducing the standard (1 -I- 3) foliation of the spacetime with 
the normals n = ~Ndt to Et where N is the lapse function. Let s be a unit normal 
to a set of time-like surfaces such that n • s = 0. We can now define two null vector 
fields by 

£ = N{n + s); k ^ {l/2N){n - s) (52) 

The primary vector field we are interested in is £ while k is an auxiliary null vector 
field with £ • k = —1 which is useful for the study of extrinsic geometry. (Of course, 
£ and k can be introduced without n and s but this allows a natural normalization.) 
We can now define a metric qab on the 2-dimensional surface St orthogonal to the n 
and s through standard relations: 

qab = gab + naUb - SaSb = gab + Lh + hka; qabi'' = = qabk'' (53) 
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The mixed tensor allows us to project quantities onto St- We can also define 
another projector orthogonal to by the definition H''^ = + k'^^b which has the 
properties 

nV'^T; ll\k' = 0; nVa=0; H^fc, = fc^. (54) 

We can now introduce the Weingarten coefficients as the projection of the covariant 
derivative Vd^"' by the definition 

x\ = u\Vdr = v^r + ibik^s/dH (55) 

which has the properties 

X\£'' = i^r-, Xabk' = 0; xabf = 0; Xabk'' = -lo^ = -PW^kb (56) 

where the surface gravity k is defined through the relation fJVjii = nii and LOa 
through the last equality. The only non-trivial result is Xabk"" — —PWjki, which 
can be proved (see e.g., Eq. (5.40) of ref. PJ]) by working out the components in 
the adapted coordinate system. Note that uja satisfies the relations uja£'^ — k and 
ujak"- — 0. We next define &ab by projecting Xmb to St- We get, on using i™'Xrnb = 
and k^Xrnb ~ ~<^f), the result: 

Qab = Q^^Xrnb = Xab + kaP^Xmb + Lk'^Xmb = Xab " ^a^b (57) 

which is essentially Eq. (|18l) expressed in the four-dimensional notation with suitable 
projection tensors. Using Eq. ([TT|) we see that 

Oab = Oba = VbL + tak'Vdb ~ h^a = C^^Vrn^n (58) 

This result shows that Qab is a natural projection of the covariant derivative Vm^n 
onto the surface St and, obviously, Qab^^ = = Qabk^ ■ The trace of Gab, denoted 
by 0, is given by 

Ql = = Va^ - n (59) 

which corresponds to Eq. (I2ip . It is also convenient to define a similar projection of 
uJa by = q^ijJa- We have 

= qluJa = -QbikmX'^a) ^ ^b ~ nkb{k^r) ^UJb + Kh (60) 

For computational purposes it is often convenient to relate Qab, etc. to more 
familiar quantities defined using the standard (1+3) decomposition of the metric. It 
can be shown that (see Section 10.2 of ref. 12 ) the following results hold which are 
often useful for explicit computation: 

eab^N{DmSn-Krnn)qaqb; 9 ^ N {D^s'' + K^ps'' ^ K) (61) 

,i = e"'\7„JnN + s''D^N-NK^fis''s>^; ^ Da\nN - Krans"'ql (62) 

These results again allow us to express the projection of Einstein's equations onto 
St- To do this one begins with the standard relation Wm^af-™' ~ ^ mP^ — Rma f™' 
and substitute for Va^™ using Eq. ([551) and for Vm£™ using Eq. ([5^ repeatedly. This 
leads, after some straightforward algebra, to the relation 

i?ma r = V„e™ + i"'W„,UJa + {k + 9)uJa ^ W a{K + 6) ~ ea„,fc"V„f " 

- (L^™fc"V„r' + V™fc" V„£™ -I- fc"V™V„r') £a (63) 
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The projection of R^aP^ using n"^ separates into two terms given by 

Rmn r'n« = -i?™„r'f" + (64) 

As mentioned in the main text, we now see that there is a projection along i itself 
and a projection to the 2-surface St- Of these, the i?„i„i'™£" will give the familiar 
Raychaudhuri equation while Rmn£™Qa lead to an equation which looks similar 
to Navicr-Stokcs equation. The Raychaudhuri equation can be obtained very easily 
by contracting Eq. ([63]) with £ and simplifying the result. We get 

Rmn ^"r = -e„„e"" + Ke- r^mO. (65) 

It is again conventional to separate out the trace of Qmn and define cr„i„ — Qmn ~ 
{l/2)qmn0 SO that we can write e™„e™" = cr™„cr™" + {1/2)6'^. 

The derivation of the Navier-Stokes like equation is more complicated. Contract- 
ing Eq. (|63l) with leads to the expression 

i?m«^™9a - q'y„^en+qy"'^r,^UJn + {K + 0)^la-Da{K + 9) 

-ea™fc"V„^™. (66) 

with (7"V„(k + 9) = Da{K + 0) where Da is the covariant derivative defined using 
the metric on St- The first and the last terms on the right hand side can also be 
combined in terms of Da using the relation: 

= qt^,Ql~Qi{k'^do+^o) (67) 

In obtaining this result we have repeatedly used the orthogonality condition 8^£" = 
= Safe" and the definitions P^Vnkm = Wm, qJ^LUm — ^j- Therefore we can combine 
the first and last terms in Eq. (155)) as 

CVmO™ - &:k"Wnirn = ^mO^ + O^f^^. (68) 

In the second term in Eq. ([55]) we introduce the Lie derivative of fin through the 
expansion 

q^e'^VmUJn = cr'v,„(i7„ - Kkn) = q:{e"'w^n^ - K^"v„fc„) 

= q2£en,,-Q^n^-Kna, (69) 

Combining all these together we get 

RmnTq: = qT£e^,n+Ona-Da{K + e)+D„,QT 

= q'X£en,n +9na-Da(^K+^^ + (70) 

where we have again set 9™ = fi™ + (l/2)(5™6'. 

This result can be re-expresscd in several different ways. For example, the con- 
vective derivative of r2„ has the form 

Dtna = q^rW^iln = qT£e^n " 6;"^!™ (71) 
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using which we can express Eq. ([70]) in terms of the convective derivative. Alter- 
natively, we can express it in the coordinate system adapted to S, in which the 
projection of the Lie derivative has the form 

q'S£e^^^^+v'' DBnA + nsDAv''. (72) 

where we have used £ = do + v^ba- Then Eq. ([70)) becomes, on using — oja for 
transverse components, because kA — 0, 

RmnrqA = (do + v"" Db)uja + DaV^ + OuJA + Dga^ -Da(^k+ ^9^ (73) 

This is the form which was used in the main text. 
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